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We discuss the combination of NNLO standard QCD evolution and scheme-invariant analysis for unpolarized DIS 
structure functions data as a method to reduce the theoretical errors on the determination of Qs(M|) to ~ 1% 
in order to match the accuracy forseen for experimental errors from future high statistics measurements. 
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1. INTRODUCTION 

The final HERA-II data on unpolarized deeply in- 
elastic scattering (DIS) structure functions, com- 
bined with the present world data, will allow to 
reduce the experimental error on the strong cou- 
pling constant, as(M|), to the level of 1% 0. 
On the theoretical side, the next-to-leading order 
(NLO) analyses have intrinsic limitations which 
allow no better than 5% accuracy in the deter- 
mination of Q!s [2]' In order to match the ex- 
pected experimental accuracy, analyses of DIS 
structure functions need then to be carried out 
at the NNLO level, which requires the knowledge 
of the /3-function and anomalous dimensions at 
the 3-loop level and the massless and massive 2- 
loop Wilson coefficients. With the recent com- 
putation of the 3-loop anomalous dimensions [2], 
the whole scheme independent set of quantities 
is known, thus allowing a complete NNLO study 
of DIS structure functions. At the same time we 
think that combining the standard QCD analysis 
and fits based on scheme-invariant evolution will 
provide a valuable tool to reduce theoretical and 
conceptual uncertainties in high-precision analy- 
ses aiming at 1% accuracy. 

Our final goal is to perform the full NNLO anal- 
ysis of DIS structure functions aiming at an high- 
accuracy determination of as and the extraction 
of a set of parton distribution functions (PDFs) 
with fully correlated errors. A complete analysis 



refers to both the singlet and non-singlet evolu- 
tion. In the present letter we will concentrate 
only on the singlet sector, referring the reader to 
the recent non-singlet analysis 0]. 

2. QCD EVOLUTION EQUATIONS 

Evolution equations of DIS structure functions 
depend, in the standard QCD approach, on two 
arbitrary scales which are introduced in the calcu- 
lation: the renormalization and the factorization 
scales. 

The renormalization scale dependence of any 
observable can be removed only summing the per- 
turbative series to all orders. Its presence is then 
unavoidable in any fixed order result. Moreover, 
the dependence of the result on the variation of 
this unphysical scale can be used to give a rough 
estimate of the theoretical error due to higher or- 
der corrections. 

On the other hand, if we consider the depen- 
dence of the result on the factorization scale we 
may follow two approaches. The first one is to 
consider the evolution of parton distribution func- 
tions which are process-independent but depend 
on the adopted factorization scheme. The sec- 
ond one is to study evolution equations for phys- 
ical observables. In these equations the role of 
anomalous dimensions for mass factorization is 
played by physical anomalous dimensions, which 
are independent of the factorization scheme but 
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depend on the process and the observables con- 
sidered. These two choices define what we call the 
standard QCD analysis and the scheme-invariant 
analysis. 

In the standard QCD analysis one introduces 
a parameterization for the different PDFs at a 
given reference scale. The PDFs are then evolved 
up to the actual scale of the process, solving the 
evolution equations for mass factorization. Struc- 
ture functions are then constructed as a convolu- 
tion of the PDFs and the corresponding Wilson 
coefficients. As a last step, a multi-parameter 
fit is then performed to extract the value of the 
coupling constant and determine the parameters 
entering the PDF parameterization. 

In a scheme- invariant analysis the parameter- 
ization of the observable at the reference scale 
Qq is extracted from the data. The value of the 
observable at the scale is determined solving 
the evolution equations with physical anomalous 
dimensions as evolution kernels. Finally a one- 
parameter fit is performed to extract the value 
of as{M|). Once the analysis is completed the 
parton densities in any factorization scheme can 
be extracted along with the respective experimen- 
tal errors. The advantage of considering factor- 
ization scheme invariant evolution equations for 
physical observables resides in the fact that the 
input distributions are observables. Full exploita- 
tion of this advantage requires therefore high 
statistics measurements to minimize errors on the 
input distributions. Furthermore, the correla- 
tions between the measured input distributions 
have to be considered in detail. 

Once more we would like to stress that the two 
analyses are complementary and not mutually ex- 
clusive. Thus performing both of them and com- 
paring the extracted values of ^^(Mf ), or corre- 
spondingly Aqcd, provides a test of stability to 
determine the QCD parameter. 



3. PHYSICAL ANOMALOUS DIMEN- 
SIONS 

When considering the singlet evolution the quark- 
singlet and gluon PDFs can be mapped into a pair 
of structure functions via the matrix of Wilson 



coefficients, C , (H] : 



(1) 



In Eq. 1^, as we will do in the following, we work 
in Mellin space, where convolutions are ordinary 
products. 

The singlet evolution equations read 



dt [ Fff 



where the evolution variable is 
/3o ""a^iQl) ■ 
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The coupling constant is related to the usual 
strong interactions coupling via the relation 
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and its expansion to 3-loops reads 
2 1 r A InL 



(5) 



where 



L = In 



^^QCD 



the (3 function is defined as 
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and, in the case of SU{2>)c, the coefficients enter- 
ing up to 3-loops are 



= ll-giV/, 

OQ 

Pi = IQ2- -Nf, 
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with Nf the number of flavors. The physical 
anomalous dimensions can be expressed in 
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terms of the anomalous dimensions and the Wil- 
son coefficients as [H] 



^ij — 
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dt 



(9) 



W{as (Q2)) 



Here denotes the unpolarized anomalous di- 
mensions which are related to the evolution ker- 
nels in X— space by 



7rrm = -2 / dx X P,nn{x) , TO, n = g, 5 (10) 

Jo 

and Cfj^ are the Mellin transforms of the Wilson 
coefRcients 



Cf„ = / dxx''-'Ci,,n{x) 
Jo 



(11) 



While the anomalous dimensions and the Wilson 
coefficients are, separately, factorization-scheme 
dependent quantities, the combinations ^ defin- 
ing the physical anomalous dimensions are fac- 
torization scheme invariants, order by order in 
perturbation theory. 

Different pairs of structure function can be 
taken into consideration, in particular: 

• F2 and dF2/dt ; 

• F2 and Fl 03 . 
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In the case of polarized DIS a combined AI S and 
scheme invariant analysis was carried out to NLO 
for the structure functions gi and dgi/dt in j7). 

Here we present the NNLO physical anomalous 
dimensions for the coupled evolution of the struc- 
ture functions F2 and dF2/dt, alongside with 
those at LO and NLO, cf. [5] : 
LO : 
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The analytic structure of the physical anomalous 
dimensions -^f"' in a;— space is difficult to obtain 
since they require the computation of the inverse 
Mellin transform of products of coefficient func- 
tions and inverse splitting functions, which are 
highly non-trivial already in the simplest cases 
(for an example see We therefore perform 

the evolution in Mellin space. The result, com- 
puted for integer N is then analytically continued 
to arbitrary complex N using representations of 
harmonic sums with arbitrary precision jJOl and 
the a;— space result is obtained through a single 
numerical contour integral. 

4. HEAVY FLAVORS CONTRIBUTION 

Heavy flavor contributions to DIS structure func- 
tions are known to be sizable in the kinematic 
region spanned by HERA. An example being the 
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Figure 1. NLO scheme invariant evolution for the 
singlet part of the structure function F2 for four 
light flavors. 
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Figure 2. NLO scheme invariant evolution for the 
singlet part of dF2/dt for four light flavors. 

6. CONCLUSIONS 



structure function F2 which, depending on actual 
event kinematics, can receive contributions from 
heavy flavors up to the level of 20 — 40%. Any 
analysis aiming at extracting as form DIS struc- 
ture functions data with an accuracy of 1% 
must, therefore, take into account heavy flavor 
contributions. Recently a parameterization of 
heavy flavor Wilson coefhcients in Mellin space 
has been derived thus allowing a direct in- 
corporation into computer codes which solve the 
evolution equation in Mellin space. 



The future high precision HERA-II data will al- 
low a reduction of the experimental error on the 
determination of as to ^ 1%. On the theoret- 
ical side, the inclusion of NNLO corrections is 
mandatory to cope with such a level of accu- 
racy. In view of a high accuracy determination of 
the strong coupling constant we think that com- 
bining the standard MS analysis and fits based 
on factorization-scheme invariant evolution could 
provide a method to have better control on theo- 
retical and conceptual errors on as. 
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5. NUMERICAL RESULTS 

While full numerical implementation of the 
NNLO scheme invariant evolution is almost fin- 
ished, as a preliminary result, in Fig. 1 and Fig. 
2 we present the scheme invariant evolution for 
the structure functions F2 and dF2/dt at NLO 
for four light flavors. In the present computation 
the initial form of the observables is built up as 
a convolution of Wilson coefficients and PDFs at 
the reference scale Qq — IGeV^, using parame- 
terization of |12|. 
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